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Abstract. The thermodynamics of an electrically charged, multicomponent fluid with spontaneous electric 
and magnetic dipoles is analysed in the presence of electromagnetic fields. Taking into account the chemical 
composition of the current densities and stress tensors leads to three types of irreversible terms: scalars, 
vectors and pseudo- vectors. The scalar terms account for chemical reactivities, the vectorial terms account 
for transport and the pseudo- vectorial terms account for relaxation. The linear phenomenological relations, 
derived from the irreversible evolution, describe notably the Lehmann and electric Lehmann effects, the 
Debye relaxation of polar molecules and the Landau-Lifshitz relaxation of the magnetisation. This formal- 
ism accounts for the thermal and electric magnetisation accumulations and magnetisation waves. It also 
predicts that a temperature gradient affects the dynamics of magnetic vortices and drives magnetisation 
waves. 

PACS. 05.70.Ln, 75.76.+j, 47.65.-d 



1 Introduction 

Spin caloritronics is mainly focused on studying the ef- 
fects of a temperature gradient on the time evolution of 
the distribution of the local spin average of a physical 
system [T]. In many experimental situations, the system 
can be treated as a classical continuum with magnetisa- 
tion on the scale of interest where the quantum fluctua- 
tions average out and the underlying microscopic struc- 
ture is smoothed out. In such as case, the local system 
is sufficiently large from a microscopic perspective to be 
treated classically, but it is sufficiently small from a macro- 
scopic perspective to be considered as infinitesimal. In 
order to understand the interplay between the magneti- 
sation and a temperature gradient, a consistent classical 
thermodynamic theory of continuous media with magneti- 
sation needs to be established. It is of interest to include 
also the electric counterpart of the magnetisation, that is 
the electric polarisation. The thermodynamics of continu- 
ous media is presented in details in textbooks by Gurtin 
et al. [?], Wilmanski [?] and Silhavy [?]. In order to ac- 
count for the intrinsic rotation of the matter, we need to 
introduce explicitly an intrinsic angular velocity. The in- 
troduction of intrinsic rotations into the thermodynamics 
of continuous media was discussed notably by Miiller [2] 
and Muschik [3,4. The electrodynamics of an electrically 
charged continuous medium interacting with clectromagt- 
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netic fields is detailled in textbooks by Eringen [?] and 
O'Dcll [?]. Some pioneering work on the non-equilibrium 
thermodynamics of a continuous medium with electro- 
magnetic fields was carried out by Liu [?] , flutter [?] and 
de Groot and Mazur [?]. In the present work, we follow 
essentially the approach of Stiickclberg detailed in ref- 
erence [5] and extend his formalism to describe locally 
the thermodynamics of an electrically charged continuous 
medium with spontaneous electric and magnetic dipoles 
in the presence of electromagnetic fields. 

In order to be able to define intensive fields such as the 
temperature, the chemical potential, the electrostatic po- 
tential, the electric field and the magnetic induction field, 
we require the local infinitesimal system to be homoge- 
neous, uniform and at equilibrium. 

We introduce phenomenological relations that express, 
in terms of the chemical constituents of the continuous 
medium, the mass density, the electric charge density, the 
intrinsic angular mass density, the electric polarisation 
and magnetisation, the momentum and the intrinsic an- 
gular momentum. Using these phenomenological relations, 
we express the electric polarisation and magnetisation cur- 
rent tensors as well as the stress and angular stress tensors 
in terms of the chemical current densities. We show that 
by expressing the state fields, such as the electric polarisa- 
tion and the magnetisation, and the dynamical fields, such 
as the momentum density and the intrinsic angular mo- 
mentum density, in terms of the different elementary con- 
stituents of the medium, three types of irreversible terms 
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appear. First, there are scalar terms that account for the 
chemical reactivities and are expressed as the product be- 
tween the chemical reaction rate densities and the chem- 
ical affinities. Second, there are vectorial terms that ac- 
count for the transport and are expressed as the dot prod- 
uct between the current densities and the forces. Third, 
there are pseudo-vectorial terms that account for the re- 
laxation and are expressed as the dot product between the 
intrinsic rotation rate densities and the intrinsic torques. 

To illustrate our formalism, we derive explicit ex- 
pressions for the irreversible relations characterising the 
Lchmann and electric Lehmann effects, and the relax- 
ation of the electric and magnetic dipoles. Furthermore, 
we show how this formalism accounts for the thermal and 
electric magnetisation accumulations and magnetisation 
waves. We also predict how a temperature gradient affects 
the dynamics of magnetic vortices and drives magnetisa- 
tion waves. 

The structure of this publication is the following. In 
Sec. [H we define the continuity equations, which describe 
the local time evolution for the material state variables. 
Then, we express the time evolution equations in terms of 
the different chemical constituents. In Sec. [31 using the fact 
that the energy is a state function, we combine the time 
evolution equations for the internal energy, the entropy 
and the other extensive material state variables. This en- 
sures the local compatibility of the first and second laws 
of thermodynamics with the laws of mechanics. Through 
symmetry arguments, we establish the thermostatic, the 
reversible and irreversible thermodynamic equations. Fi- 
nally, in Sec. 21 we establish the irreversible thermody- 
namic phenomenology and describe some physical phe- 
nomena to illustrate our formalism, with a particular em- 
phasis on spintronics. 



2 Thermodynamic description 

2.1 Continuity equations 

The state of a continuous medium is defined by a set of 
matter state variables. The local state of a continuous 
medium is defined by a set of matter state fields that are 
function of the space and time coordinates. To keep the 
notation concise, we do not explicitly write the space and 
time dependence of the state fields. The local thermody- 
namic state of a classical electrically charged continuous 
medium consisting of N chemical substances with electric 
polarisation and magnetisation interacting with electro- 
magnetic fields is determined by the following state fields: 

• the entropy density field s, 

• the number density fields ua of ./V chemical substances, 
where A 6 {1,...,N}, 

• the electric charge density field q, 

• the spontaneous electric polarisation field P, 

• the spontaneous magnetisation field M, 

• the velocity field v, 

• the intrinsic angular velocity field u). 



It is important to stress that the system does not include 
the electromagnetic fields in the vacuum, but only their 
interaction with the matter fields. Thus, the state fields q, 
P and M, accounting for the electromagnetic properties 
of the matter, are purely matter fields, in contrast to our 
previous work where the electric displacement field 
D and the magnetic induction field B were accounting for 
the electromagnetic properties of the matter and of the 
electromagnetic fields in the vacuum. 

Moreover, in this publication, we restrict our analysis 
to a system with spontaneous electric polarisation P and 
magnetisation M, whereas in our former work we 
considered only a system with induced electric polarisa- 
tion and magnetisation. 

The classical time evolution of every extensive phys- 
ical property in a local infinitesimal system is expressed 
by a continuity equation, which is a local detailed balance 
equation. In this publication, we shall refer to the conti- 
nuity equations expressed in terms of the material time 
derivative as the material continuity equations. 

The material continuity equation for an extensive 
scalar property F is given by [5J, 

/ + (V-v)/ + V-j / =p / , (1) 

where / is the material time derivative of the scalar den- 
sity f, jf is the vectorial diffusive current density and pf 
is the scalar source density of F. Note that the frame- 
independent material time derivative operator is related 
to the partial time derivative operator by, 

• = d t + (v • V) . 

The material continuity equation for an extensive vec- 
torial property F is given by [5J, 

f +(V-v)f + V-j f = p f , (2) 

where f is the material time derivative of the vectorial 
density f , jf is the rank-2 tensorial diffusive current density 
and pf is the vectorial source density of F. 

The material continuity equation ([1} for the entropy 
yields, 

s + (V-v)s + V-j s =p s >0 , (3) 

where s is the entropy density, j s is the diffusive entropy 
current density and p s is the entropy source density that is 
positive-definite in order to satisfy locally the second law 
of thermodynamics. Note that the strictly positive part, 
i.e. p s > 0, accounts for the irreversibility by breaking the 
symmetry of the dynamics under time inversion, which 
defines a time arrow. 

The material continuity equation ([l} for the number 
of elementary units of chemical substance A is given by, 

hA + ( V ■ v) n A + V ■ j A = ^ VaA UJa ' ( 4 ) 

a 

where n A is the number density of the chemical substance 
A, j A is the diffusive current density and the source density 
Pa is expressed explicitly as the sum over all the chemical 
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reactions a of the product of the stoichiometric coefficients 
v a A and the reaction rate densities u a . 

The material continuity equation (JT|) for the electric 
charge is given by, 

q+(V-v)q + V-j, = , (5) 

where q is the electric charge density, j 9 is the diffusive 
electric current density and there is no electric charge 
source density p q since the electric charge is an invari- 
ant. Note that the electric current density j is the sum of 
the convective and diffusive electric current densities, i.e. 

j = gv +j g . (6) 

The material continuity equation (|2|) for the sponta- 
neous electric dipoles is given by, 

P + (V-v)P + V-j P = ^ flJfxp A , (7) 

A 

where P is the spontaneous electric polarisation, jp is the 
diffusive electric polarisation current tensor. The source 
density pp of the spontaneous electric dipoles of the 
chemical substances A is an irreversible term accounting 
for the relaxation of the spontaneous electric dipoles pa- 
It is due to the rotational motion of the chemical sub- 
stances A, that is expressed in terms of the intrinsic rota- 
tion rate densities fl ™ of the matter where the superscript 
m stands for "matter" . 

The uniformity of the local infinitesimal system implies 
that there is a unique local rotation rate density for 
each chemical substance .4. that is defined with respect to 
the local frame where the electric polarisation P does not 
rotate. This frame corresponds to the local frame where 
the intrinsic rotation of the matter vanishes, i.e. uj = 0. 
since the spontaneous electric polarisation is rigidly at- 
tached to the matter. 

The material continuity equation ^ for the sponta- 
neous magnetic dipoles yields, 

M + ( V • v) M + V • j M = 

[lA n A (m, xB)|f2fx m A ) , ( 8 ) 

A 

where M is the spontaneous magnetisation, jM is the dif- 
fusive magnetisation current tensor. The source density 
Pm of the spontaneous magnetic dipoles of the chem- 
ical substances A consists of two terms. The first term is a 
reversible term accounting for the precession of the mag- 
netic dipoles around the magnetic induction field B 
that is expressed in terms of their gyromagnetic ratios 7,4- 
The second term is an irreversible term accounting for the 
relaxation of the magnetic dipoles with respect to the 
magnetic induction field B that is expressed in terms of 
the intrinsic rotation rate densities fi^ 1 where the super- 
script M stands for "magnetisation" . These rotation rate 
densities consist of two contributions. The first contribu- 
tion corresponds to the micromagnetic hypothesis [7] , that 
accounts for the fact that electrons are point-like particles 



that have no known substructure and thus no intrinsic 
angular mass. The second contribution accounts for the 
rotational motion of the chemical substances A. The sec- 
ond contribution is generally much smaller than the first. 

The uniformity of the local infinitesimal system implies 
that there is a unique local rotation rate density fl^ 1 for 
each chemical substance A, that is defined with respect 
to the local frame where the magnetisation M docs not 
rotate. In contrast to the case of the electric polarisation, 
this frame docs not correspond to the local frame where 
the intrinsic rotation of the matter vanishes, since the 
spontaneous magnetisation is not rigidly attached to the 
matter. Note that by symmetry, only axial vectors such 
as the magnetic dipoles are allowed to precess. Thus, 
no precession of the electric dipoles pa is to be expected. 

The material continuity equation for the momen- 
tum is given by, 

p+(V-v)p- V<x = f cxt , (9) 

where p is the momentum density, cr is the stress tensor 
and f oxt is the external force density. The momentum den- 
sity p {tia, v) of the matter is proportional to the velocity 
v, i.e. 

p(nA,v) = m(n A ) v , (10) 

where m (tia) is the mass density. 

The material continuity equation ([T]) for the mass 
yields, 

m+(V-v)m = 0. (II) 

There is no diffusive mass current density j m by definition 
of the local centre of mass and there is no mass source 
density p m since the mass is a non-relativistic invariant. 
Substituting the relation pTJ|) and the continuity equa- 
tion (jTTJ) for the mass into the continuity equation © for 
the momentum, the latter yields Newton's second law of 
motion, i.e. 

77iv = f cxt + V-cr . (12) 

The material continuity equation ^ for the intrinsic 
angular momentum yields, 

s + (V-v)s- V-@ = r oxt , (13) 

where s is the intrinsic angular momentum density, is 
the intrinsic angular stress tensor and r cxt is the intrinsic 
torque density. The local infinitesimal system is consid- 
ered as a homogeneous sphere. Thus, the intrinsic angular 
momentum density s (tia, of the matter is proportional 
to the intrinsic angular velocity <*>, i.e. 

s(ru,w) = I {tia) u , (14) 

where / {tia) is the intrinsic angular mass density, which 
is a scalar density. For a homogeneous sphere of constant 
infinitesimal radius dr, the intrinsic angular mass density 
/ is related to the mass density m by, 

I = - m dr 2 . (15) 
5 
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Thus, the material continuity equation (jllj) for the 
mass and the relation (|15|) imply that the material conti- 
nuity equation ([1} for the intrinsic angular mass is given 
by 

7+(V-v)/ = 0. (16) 

There is no diffusive intrinsic angular mass current density 
ji and no intrinsic angular mass source density pj since 
for a homogeneous local system the intrinsic angular mass 
density / is totally determined by the mass density m. 

Substituting the relation (Ti"4")) and the continuity equa- 
tion (|16p for the intrinsic angular mass into the continuity 
equation p3[) for the intrinsic angular momentum, the lat- 
ter yields Newton's second law in intrinsic rotation, i.e. 

Iu> = T CKt + V . (17) 

The material continuity equation (|T|) for the energy is 
given by, 

e+(V-v)e + V-j e = f cxt -v + T oxt -u; , (18) 

where e is the energy density, j e is the diffusive energy 
current density and the energy source density p e is the 
power density due to the external force densities f cxt and 
external intrinsic torque densities T ext . The energy den- 
sity e (s, tia, q, P, M, v, u>) is the sum of the translational 
kinetic energy density, i.e. i m v 2 , the intrinsic rotational 
kinetic energy density, i.e. 1 I uj 2 , and the internal energy 
density, i.e. u (s, ha, q, P, M), 

e(s,n A ,q,P,M,v,u) = ^mv 2 + ^Iv 2 

+ u(s,n A ,q, P,M) , 

where the internal energy density is the energy density in 
the local rest frame where v = and u> = 0. Using the 
continuity equations for the mass (jllj) and the intrinsic 
angular mass (| , Newton's second law in translation (|12p 
and in intrinsic rotation (|17[) . the time derivative of the 
relation (fl"9| yields, 

e = u + (V • v) (u - e) + (V • a) ■ v + (V • 0) ■ u> 
+ f oxt • v + x cxt • w . (20) 

Substituting the relation (|20j) into the continuity equa- 
tion (fT5|) for the energy and using the vectorial identities 

(V • er) • v = V CT ■ (er • v) - er • (V v) , 
(V • 0) ■ u> = V© • (0 ■ u) - ■ (V © u) , 

yields the material continuity equation ([T} for the internal 
energy density, i.e. 

u+(V • v)it + V-j u = <7-(V©v) + <9-(V0u;) , (21) 

where the symbol denotes a symmetrised tensorial prod- 
uct and the indices CT or @ denote that there is a dot 
product between the covariant differential operator V and 
the contravariant components of the stress tensors er and 
respectively. 



The internal energy current density j u is related to the 
energy current density j e by, 

j e =j«-crv-6»w. (22) 

Note that the source density terms cr • ( V v) and 
■ (V uj) in the continuity equation (|21[) for the inter- 
nal energy account for the elastic expansion or contraction 
of the matter and for the dissipative translational and ro- 
tational viscous friction. The dissipation is defined as the 
irreversible part of the internal energy source density p u . 



2.2 Dynamics in terms of the chemical composition 

2.2.1 Time evolution of the electric charge density 

The electric charge density q is defined as the product of 
the number density n A and the electric charge q A carried 
by the elementary units of chemical substances A [6], i.e. 

q = 2J n A q A . (23) 
A 

In order to characterise physically the time evolution of 
the electric polarisation, we substitute the expression ([23)) 
for the chemical composition of the electric charge density 
and the continuity equation for the chemical substance ^ 
into the continuity equation §5§ for the electric charge. The 
latter then becomes, 

X! (- ( V '.U) + "a A u a q A ) + V-j, = . (24) 

A. a 

Using the fact that the electric charge q A is an invariant, 
i.e. 

(V • j A ) q A = V • (q A j A ) , 

the continuity equation |24|) for the electric charge can be 
recast as, 

X v aA w a q A + V • (j q - X q A ja) = . (25) 

a, A A 

The continuity equation (|2"5]) has to hold for any current 
flow, which yields an explicit expression for the diffusive 
electric current density, i.e. 

h = ^2 IaU ■ (26) 

A 

Moreover, it has to hold for every chemical reaction a, 
which implies that 

V «A q A = , (27) 

A 

and means that the chemical reaction a preserves the total 
electric charge in the local infinitesimal system. 
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2.2.2 Time evolution of the electric dipoles 



2.2.3 Time evolution of the magnetic dipoles 



The electric polarisation P is defined as the product of the 
number density ha and the electric dipoles pa carried by 
the elementary units of chemical substances A, i.e. 



P = ^2 n A p A 

A 



(28) 



In order to characterise physically the dynamics of the 
electric polarisation, we substitute the expression (|2"5)l for 
the chemical composition of the electric polarisation and 
the continuity equation for the chemical substance $3) into 
the continuity equation ([7]) for the electric dipoles. The 
latter then becomes, 



E (™ Pa - (V ' 3a) Pa + ^2 VaA u a Pa) + V • ji 



V £2% x p A 



A 



(29) 



Using the vectorial identity, 

(V -u) Pa = Vj ■ (p A Ja) ~ (U ■ V)pa , 

where the index j denotes that there is a dot product be- 
tween the covariant differential operator V and the con- 
travariant current density j a, the continuity equation (|29[) 
for the electric dipoles can be recast as, 



E 

A 



(n A P A - ft a x Pa + (jA ■ V) p A 



^ V a A w« PA + Vj ■ (jp - ^ PA JA ) 

a, A 



A 



(30) 



The continuity equation (|30|) has to hold for any current 
flow, which yields an explicit expression for the electric 
polarisation current tensor, i.e. 



Jp = E PA J y 



(31) 



The condition ([3"T]) implies that the continuity equa- 
tion (l30l) reduces to, 



E ( U A PA- ft A XPA + (jA • V)p A + ^2 VaAUJaPA^ = , 

A a 

. (32) 

which describes the time evolution of the electric dipoles 
Pa of all the chemical substances. 

The local time evolution of the electric dipoles pa of 
a specific chemical substance A is given by, 

PA = -PA* ft™- (JA-V)PA- VaAUaPA , (33) 



"A 



where the first, second and third terms on the RHS de- 
scribe respectively the relaxation, the transport and the 
chemistry of the electric dipoles pa- Note the intrinsic 
rotation rate densities £2™ are functions of the electric 
dipoles of all the chemical substances. They account for 
the dissipative couplings between the electric dipoles of 
the different chemical substances. 



The magnetisation M is defined as the product of the 
number density ua and the magnetic dipoles mj carried 
by the elementary units of chemical substances A, i.e. 



M = n A m A ■ 
A 



(34) 



In order to characterise physically the time evolution of 
the magnetisation, we substitute the expression (f3"4"| for 
the chemical composition of the magnetisation and the 
continuity equation for the chemical substance ((4|) into 
the continuity equation 1(5)) for the magnetic dipoles. The 
latter then becomes, 

Y2 ( n A «ia - ( V • j A ) hia + ^2 VaA m A ^j + V ■ j M 

A a 

= ^2 wA n A m A x B + £2™ x m A J • (35) 



Using the vectorial identity, 

(V • ja) m A = Vj ■ (mA jA) - (jA ■ V) m A , 

the continuity equation (|35[) for the magnetic dipoles can 
be recast as. 



E 

a, A 



The continuity equation (J36J) has to hold for any current 
flow, which yields an explicit expression for the magneti- 
sation current tensor, i.e. 



E { UA ( m -4 - 7A mA x B) - £2™ xm A + (jA • V) m A 

A 

VaA u a m A + Vj • (j M - ^ m A Ja) = . 

(36) 



JM = 2^ 

A 



mA Ja 



(37) 



^ (n A (riiA - 7a mA x B) - £2™ x m A 

A 

+ (ja • V) m A + ^2 VaA u a mA J = 



Note that this expression for the magnetisation current 
tensor is the classical counterpart of the tensorial "spin- 
current" used in spintronics [5] . The condition (|37[) implies 
that the continuity equation (|36p reduces to, 



(38) 



which describes the local time evolution of the magnetic 
dipoles of all the chemical substances. 

In order to establish the local time evolution of the 
magnetic dipoles mA of a specific chemical substance A 
in the presence of other chemical substances B carrying 
magnetic dipoles mg. wc need to take explicitly into ac- 
count the local magnetic torque that couples mA and m^. 
This torque is antisymmetric under the permutation of dif- 
ferent chemical substances A and B, which implies that 
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the torque density vanishes after summation over all the 
chemical substances, i.e. 



^ Iab n A n B (m A x m B ) = , 



(39) 



A,B 



where "f AB & the symmetric coupling coefficient. Note that 
the magnetic torques do not affects the magnetisation M, 
which is denned as the density of magnetic dipoles of all 
the chemical substances. Thus, they do not affect the Lar- 
mor energy of the local system. Note that by symmetry, 
the interaction torques are allowed only for axial vectors 
such as the magnetic dipoles and mj. Thus, no in- 
teraction torque of the electric dipoles and p B is to be 
expected. 

The equations (f38|) and (|39"1) imply that the local time 
evolution of the magnetic dipoles of a specific chemical 
substance A is given by, 

n A rh A = 7a n A (m A x B) - x fl™ - (j A ■ V) m A 



v aA u a m A + ^ Iab n A (m A x n B m B ) 



(40) 



where the first, second, third and fourth terms on the RHS 
describe respectively the precession, the relaxation, the 
transport and the chemistry of the magnetic dipoles m^, 
and the last term describes the local interaction of the 
magnetic dipoles with the magnetic dipoles m B of 
the other chemical substances B. Note the intrinsic ro- 
tation rate densities fl™ are functions of the magnetic 
dipoles of all the chemical substances. They account for 
the dissipative couplings between the magnetic dipoles of 
the different chemical substances. 



2.2.4 Time evolution of the mass density 

The mass density m is defined as the product of the num- 
ber density n A and the mass m A carried by the elementary 
units of chemical substances A [5], i.e. 



m = n A m A 
A 



(41) 



In order to characterise physically the time evolution of 
the mass density, we substitute the expression (j4i~j) for the 
chemical composition of the mass density and the continu- 
ity equation for the chemical substance (j3]) into the conti- 
nuity equation for the mass. The latter then becomes, 

X] ( - ■ Ja) m A + ^2 Va A U>am A \ = 



(42) 

Using the fact that the mass m A is a non-relativistic in- 
variant, i.e. 

(V • ja) m A = V ■ {m A j A ) , 
the continuity equation (|42|) for the mass can be recast as, 

^ VaA QJg m A - V • ( m A j A ) = . (43) 



The continuity equation (|43|) has to hold for any current 
flow, which defines a set of frames where the local centre 
of mass of the matter element is at rest, i.e. 



^2 m A j A = 



(44) 



Moreover, it has to hold for every chemical reaction a, 
which implies that 



^ v aA m A = 



(45) 



and means that the chemical reaction a preserves the total 
mass in the local infinitesimal system. 



2.2.5 Time evolution of the momentum 

The momentum density p = m v is defined as the density 
of momentum carried by the elementary units of chemical 
substances A of velocities v A , i.e. 



P = ^ n A m A v A 

A 



(46) 



It consists of convective part and a diffusive parts accord- 
ing to, 

2J n A m A v A = 2J n A m A v + ^ m A '$ A . (47) 

A A A 

Thus, using the relation (|4"4l . the velocity of the matter 
element is found to be, 

v = n A mj^j ^2 n A ™>A ^A^j , (48) 

and represents the velocity of the local centre of mass. 

In order to characterise physically the time evolution 
of the momentum, we substitute the expression (|4*6")l for 
the chemical composition of the momentum density and 
the continuity equation for the chemical substance ((3]) into 
the continuity equation ([§]) for the momentum. The latter 
then becomes. 



^ [n A m A w A - (V ■ j A )m A v A + ^ v aA w a m A x A ^ 

A a 

- V • a = f cxt . (49) 
Using the vectorial identity, 

(V • j A ) v A = Vj • (v A Qj A ) ~ (j A -V)v A , 
and the extensivity of the force, i.e. 

f ext = f% xt , 
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where f| represents the external force density acting on 
the chemical substance A, the continuity equation (|49| for 
the momentum can be recast as, 

(n A m A v A + m A (j A ■ V) v A - f| xt ) (50) 

A 

+ E VaA UJa mA - Vj • (cr + ^ m A v A j A ^j = . 

a, A A 

The matter stress tensor er is split into a reversible part 
due to the pressure P and an irreversible part cr according 
to, 

cr = -P1 + & , (51) 

which implies that the continuity equation (|50p is recast 
as, 

5^ (n A m A v A + m A {u ■ V) v A - fT*) + V P (52) 
+ ^ v aA u: a m A v A - Vj • (a- + ^ to^ j A J = . 

a, A A 

The continuity equation (|52"j) has to hold for any cur- 
rent flow, which yields an explicit expression for the irre- 
versible part of the stress tensor, i.e. 



^ m A v A j A 

A 



(53) 



Moreover, every chemical reaction preserves the total mo- 
mentum in the local infinitesimal system, i.e. 



^ v aA m A v A = 



(54) 



Taking into account the fact that the pressure P of the 
continuous medium is the sum of the partial pressures P A 
of the different chemical substances A multiplied by the 
dimensionless fugacity coefficients (p A [?], i.e. 



P = J2 <PA P A , 



(55) 



and using the conditions ([53]) and ((54)) . the continuity 
equation (|5"2")) reduces to, 

(^A m A v A +m A (j A ■ V) Vj4 - f% xt +V (tp A P A ) \ = . 

A 

(56) 

which describes the local time evolution of all the chemi- 
cal substances. Note that the dimensionless fugacity coef- 
ficient tp A in the relation (|55p accounts for the local inter- 
actions between the constituents of the chemical substance 
A. 

In order to establish the local time evolution of a spe- 
cific chemical substance A in the presence of other chem- 
ical substances B, we need to take explicitly into account 
the local internal force densities ig%. A exerted by the 
chemical substances B on the chemical substance A. These 
force densities do not change the local internal energy 



of the system, which means that they are the densities 
of conservative forces accounting for elastic collisions or 
scattering. Newton's third law implies that the sum of the 
internal force densities over all the chemical substances 
vanishes, i.e. 



E 

A,B 



f int 
I B->A 







(57) 



The equations (|56")l and ([57)) imply that the local time 
evolution of a specific chemical substance A is given by, 



n A m A v A 



fl xt - V(^P A )-m A (j A -V)v A + , T f^ nt 



B 

(58) 

where the first, second and third terms on the RHS de- 
scribe respectively the action of the external forces, of the 
fugacity gradient and of the transport, and the last term 
describes the action of the internal forces due to the lo- 
cal interaction of the chemical substance A with the other 
chemical substances B. 



2.2.6 Time evolution of the intrinsic angular momentum 

The intrinsic angular momentum density s = I uj is de- 
fined as the density of intrinsic angular momentum car- 
ried by the elementary units of chemical substances A of 
intrinsic angular velocities u} A , i.e. 



s = 2J n A I A lj a , 

A 



(59) 



where n A I A represent the intrinsic angular mass density 
of the elementary units of chemical substance A and the 
relations (jT5J) and (|4"T)) imply that 



n A I A — ~ n A m A dr 2 
5 



(60) 



The intrinsic angular momentum density s consists of con- 
vective and diffusive parts according to, 



} J n A I A oj a = J~] n A I A u> + ^ I A Q 



A J 



(61) 



where the sum of the diffusive intrinsic angular momenta 
of all the chemical substances A vanish, i.e. 



(62) 



which defines a set of frames where the matter has no aver- 
age intrinsic rotational motion. Thus, the intrinsic angular 
velocity of the matter element is found to be, 



n A I A ^j ^2 n A I A Ujij , 



(63) 



and represents the average intrinsic angular velocity 
around the local centre of mass. 
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In order to characterise physically the time evolution 
of the intrinsic angular momentum, we substitute the ex- 
pression (|59p for the chemical composition of the intrinsic 
angular momentum density and the continuity equation 
for the chemical substance @ into the continuity equa- 
tion (|13[) for the intrinsic angular momentum. The latter 
then becomes, 



collisions or scattering. Newton's third law in intrinsic ro- 
tation implies that the sum of the internal torque densities 
over all the chemical substances vanishes, i.e. 



A.B 



int 



(69) 



^ [riA I A WA- (V • j A ) Ia^A+'^ / VaA Wq I A ^ A 



- V = T CXt . 

Using the vectorial identity, 

(V ■ ja) u a = Vj • (ui Ja) " CU ■ V) Ua , 
and the extensivity of the torque, i.e. 

ext \ _ ext 

= T A > 



(64) 



The equations (|6"5)l and (|6"9")l imply that the local intrinsic 
time evolution of a specific chemical substance A is given 

by 



n A I A U A 



I A (}A ■ V) UJ A 



B 



int 
r B->A ' 



(70) 



where the first and second terms on the RHS describe 
respectively the action of the external torques and of the 
transport, and the last term describes the action of the in- 
ternal torques due to the local interaction of the chemical 
substance A with the other chemical substances B. 



where r| xt 



represents the external torque density acting 
on the chemical substance A, the continuity equation (|64p 
for the intrinsic angular momentum is recast as, 



^] (nA I A u A + I a 0a -V)u> a - t 



ext 
A 



VaA U a Ia^ A - Vj • [0 + I A UJ A Qj 



(65) 



a, A 

The continuity equation (|65|) has to hold for any cur- 
rent flow, which yields an explicit expression for the irre- 
versible part of the intrinsic angular stress tensor, i.e. 



= - ^2 IauaOJa 



(66) 



3 Thermostatics and Thermodynamics 



3.1 Thermostatic equation, reversible and irreversible 
thermodynamic equations 



The thermostatics and thermodynamics of the contin- 
uous medium are contained within the internal en- 
ergy balance (|2"Tj) . Since the internal energy density 
u (s, n A , 5, P, M) is a state function in the local rest frame, 
the time derivative of the internal energy density field is 
given by, i.e. 



Moreover, every chemical reaction preserves the total in- 
trinsic angular momentum in the local infinitesimal sys- 
tem, i.e. 

^2 V aA I A uj A = . (67) 

A 

The conditions and (|6"T)) imply that the continuity 
equation (|65p reduces to, 

J2 (n A I a "A + I a (U ■ V) u> A - r A cxt ) = , (68) 



which describes the local intrinsic time evolution of all the 
chemical substances. 

In order to establish the local intrinsic time evolution 
of a specific chemical substance A in the presence of other 
chemical substances B, we need to take explicitly into ac- 



Ts 



B—>A 



count the local internal intrinsic torque densities r 
exerted by the chemical substances B on the chemical sub- 
stance A. These torque densities do not change the local 
internal energy of the system, which means that they are 
the densities of conservative torques accounting for elastic 



E 

A 



[lATlA + Vq - E P - B-M 



(71) 



where the temperature T, the chemical potential \ia of 
the substance A, the electric potential V, the opposite of 
the electric field — E and the opposite of the magnetic in- 
duction field — B are defined as the intensive conjugate 
fields of the extensive state fields s, tia, q, P and M re- 
spectively [9"lfTU]. i.e. 



ds 



l±A 



du 
dn A 



V = — 

dq 



du 
dP 



E = — , B = 



du 
dM 



(72) 



Note that these intensive fields can only be defined locally 
if the system satisfies the local equilibrium hypothesis. 

Using the relation (|62[) . the continuity equations for 
the entropy ([3j, the density of the chemical substance 
A ((3]), the electric charge density (O, the electric polarisa- 
tion Q and the magnetisation @ , the continuity equation 
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for the internal energy (p?Tj) is recast as, 

T(p„- (V-v)s- V-j.) 
+ ^ [i A ( ^2 VaA u a - ( V ■ v) n A - V • Ja 



polarisation and magnetisation, we deduce the identities, 



A a 

+ V( - (V-v)q- V-j, 

- E n?xp A - (V-v)P- V-jp) (73) 

A 

- B • ( lA n A (m A x B) + fif x mx) 

A 

- (V-v)M- V-j M 
+ (V • v) u + V ■ j„ = <r • (V © v) + ■ (V w) 

Using the vectorial identities, 

T(V-j s ) = V-(Tj s )- j s - vr, 
VA (V • Ja) = V • Qi A j A ) - U ■ V/m , 
V{V-j q ) = V-(j q V)-j q -VV, 
E • (V • jp) = V • ( JP • E) - j P • (V E) , 
(V • j M ) = V • Om • B) - j M • (V B) , 
(fljj 1 x PA ) = f2^ ■ (pa x E) , 
(7A n A (m A x B) ) = , 
(f2f x m A ) = C2f ■ (m A x B) , 



q V = IA V n A , 

P • E = J2 (Pa ■ E) n A , 

M B = ^2 ( m A ■ B) n A ■ 



(76) 



Moreover, using the expressions (f26")l . (f5T|) and (|3T|) for the 
diffusive electric current density vector, the electric polari- 
sation and the magnetisation current tensors respectively, 
we obtain the identities, 



ii v = qA V ^ A ' 

A 

j P • E = J2 (PA ■ E) Ja , 

A 

jM • B = ^2 ( m A ■ B) Ja 



(77) 



the splitting (|5f [) of the stress tensor into a reversible 
and an irreversible part, and expressing the chemical irre- 
versibility in terms of the chemical affinities A a , i.e. 

A a = ^2 VaA (~Va) , (74) 

A 

the internal energy balance equation (|73p is recast as, 
(^u-Ts + P -J2^ A n A - qV + P - E + M-B^j (V-v) 

+ V- Tj s - Yj^ASA- jg^+jp-E+jM-B^ 

+ T Ps -^2oj a A a +j s 'VT + ^jA ■ - jq VV 

a A 

- ct-(V©v)- ©-(V©w) (75) 

- j P • (V0E) - j M • (V0B) 

- J2 ■ (PA x E) - J2 °A ■ (mA x B) = . 



Furthermore, using the expressions (p)3")) . (|rJrJ|) . (|3"Tj) 
and (|37p . for the momentum and intrinsic angular mo- 
mentum stress tensors as well as for the electric polari- 
sation and the magnetisation current tensors respectively, 
we obtain the identities, 



<T ■ 


(V0v) = 


-^2iA- 






A 


& 


(V©w) = 








A 


jP 


(V0E) = 


= Y2* a 






A 


JM 


• (V B) 


= Y,i A 



m A va 



Vv) , 
I a u) A Vuj) , 



(78) 



At this point it is useful to introduce the generalised chem- 
ical potential p, A that is defined as, 



/iA = Ma + qAV - pa E m A ■ B 



(79) 



Using the expressions (j!?3"|) . and (|34[) for the chemi- 
cal composition of the electric charge density, the electric 



Using the identities f76] ) -(|78 p and the definition f79j), 
the internal energy balance equation (|75[) is recast explic- 
itly in terms of the physical properties of the chemical 
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components, i.e. 



U- Ts + P-^ VAn^J (V-v) 



T Ps -J2^aA a - j s ■(-VT) 

a 

J2 u • ( - v^a -?aV^ 

^ Ja • ( - rru v A V v - I A u A V a; 

A 

X! J.4 ' (pa V E + m A V B 

A 

^ tf£ • (PA x E) - J2 n A ■ (mA x B) = 



(80) 



Firstly, the internal energy balance equation (|80|) has 
to hold locally for all flows. This implies that the terms 
in the first brackets have to vanish, which yields the ther- 
mostatic equilibrium equation for the continuous medium, 
i.e. 

u = Ts-P+J2 [pA+qAV- p A -E- m A -B^ n A . (81) 

A 

Secondly, the internal energy balance equation (|80|) has 
to hold locally for all currents. This implies that the terms 
in the second brackets have to vanish, which yields the 
reversible thermodynamic evolution equation for the con- 
tinuous medium, i.e. 

U = T j s + (ma + 9A V - Pa • E - m A • b) j A . (82) 

A 

Thirdly, the thermostatic equilibrium equation (|81D 
and the reversible thermodynamic evolution equation (|82[) 
imply that the internal energy balance equation (|80[) yields 
the irreversible thermodynamic evolution equation for the 
continuous medium, i.e. 

I a 

A 

+ ^ JA • ( - m A va V v - I A u) A V u^j 

A 

+ Y1 3A ■ (pa VE + m A VB) (83) 

A 

+ X^-(PAxE)+X^.(mAxB)l . 

A A J 

Finally, note that the time derivative of the thermo- 
static equation (jBTj) and the expression ([7T1) for the time 



derivative of the internal energy determines a generalised 
Gibbs-Duhcm relation that yields the time evolution of 
the intensive fields (|72p . i.e. 



sT - P 



A 



n A (jL A + < ? V r -P-E-M-B = 0. (84) 



4 Thermodynamical phenomenology 

4.1 Linear phenomenological relations and Onsager 
matrix elements 

In order to deduce the linear phenomenological relations, 
we extend the approach developed by Onsager [TIHT!?] 
to include intrinsic rotations. The expression ([53]) can 
be formally split into a scalar, a vectorial and a pseudo- 
vectorial part, which arc irreducible representations of the 
Euclidean group and have different symmetries. Thus, the 
entropy source density (|83[) is expressed formally as, 



1 

Ps = j, 



JaA a + Y,i<*-F a +Y, ^a'T; 



,1. i 



(85) 

where u> a is a scalar chemical reaction rate density, A a is 
a scalar chemical affinity, j Q is a vectorial current density, 
F a is a vectorial force, fl\ is a pscudo-vectorial intrinsic 
rotation rate density and T A is a pseudo- vectorial intrinsic 
torque. 

In the relation (|55|) . there are two types of vectorial 
current densities and forces (i.e. a <E { s ,a})- First, there 
is a thermal current density j s and an irreversible thermal 
force F s = — V T. Second, there are current densities ja 
for the chemical substances A and irreversible forces Fa 
acting on these substances. The expression for the forces 
Fa is given by the irreversible relation (|53")l . i.e. 



Fa = — V/iA — <3A V V — m A w A V v — I A oj a V u) 
+ pa V E + m A V B , (86) 

where the first term on the RHS is the chemical force, the 
second term is the Coulomb force, the third and the fourth 
terms are the viscous forces, the fifth term is the electric 
polarisation force and the sixth term is the magnetisation 
force. Note that in a stationary regime, the viscous forces, 
the electric polarisation force and the magnetisation force 
generate dielectrophoresis and magnetophoresis. 

In the relation ([55)1 , there are also two types of pseudo- 
vectorial intrinsic rotation rate densities and torques (i.e. 
A e {a'A 1 }) associated to the chemical substances A. 
First, there are intrinsic rotation rate densities fl A and 
irreversible torques T™ associated with the intrinsic ro- 
tational motion of the matter. Second, there are intrinsic 
rotation rate densities fi™ and irreversible torques 
associated with micro-magnetism and matter. The expres- 
sion for the torques T A is given by the irreversible rela- 
tion d83j), i.e. 



T™ = Pa x E , 
= mA x B 



(87) 
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The local expression of the second law requires 
the entropy source density to be locally positive defi- 
nite, i.e. p s 0. Extending Onsager's approach, in the 
neighbourhood of a local thermodynamic state where the 
scalar affinities A a , the vectorial forces F a and the pseudo- 
vectorial intrinsic torques are sufficiently small, the 
entropy source density can be expressed formally as the 
sum of quadratic forms of A a , F Q and T^, which ensures 
that it is positive definite, i.e. 



= t ( LabAa Ab + H la ? ' ( Fq F ^ 

b a, P 

£ l V(tAot£))^o, 



a, b 



A, B 
i,j 

where the phenomenological components are the Onsager 
matrix elements, which are of two different types: scalars 
L a p and rank-2 tensors \- a p and \- AB - The symmetries of 
the Onsager matrices are given by the Onsager reciprocity 
relations, i.e. 



L ab (s, n A ,q, P,M) = 
U/3 (s, n Al q, P,M) = 

L% B (s,n A ,q,P,M) 



L ba (s,n A ,q,P, -M) , 
Lp a {s,n A ,q,T>,-M) , (89) 
= \- AB (s,n A ,q,P,-M) , 



which cannot be derived within a thermodynamic ap- 
proach but require a statistical treatment since they are 
a consequence of the time reversibility of the microscopic 
dynamics [11]. The inequality ([88)) has to hold for each 
part, which implies that each quadratic form has to be 
positive definite. 

Thus, the chemical reaction rate densities u> a are re- 
lated to the chemical affinities Ab through scalar linear 
phenomenological relations, i.e. 



where the Onsager matrix has to satisfy, 



1 



(90) 



(91) 



y {Lab} > , 



in order for the scalar quadratic form in the relation 
to be positive definite. The scalar linear phenomenolog- 
ical relations ((90)) account for the irreversibility due to 
the chemistry. Note that these linear relations are only 
accurate provided that the chemical affinities Ab are suf- 
ficiently small compared to the thermal excitation, i.e. 
Ab <C ksT where ks is Boltzmann's constant, as pointed 
out by Glansdorf and Prigogine [?] and de Groot and 
Mazur [?]. 

Similarly, the vectorial current densities jo, are related 
to the vectorial forces F^ through vectorial linear phe- 
nomenological relations, i.e. 



Ja = 



(92) 



where the Onsager matrix has to satisfy, 

>0 , 



(93) 



in order for the vectorial quadratic form in the rela- 
tion ([88)) to be positive definite. The vectorial linear phe- 
nomenological relations (|9"2")l are expressed explicitly in 
terms of the currents densities j s and ] A and forces F s 
and F A as, 



'j s = L ss .(-VT) + ^L sB .F B 

B 

U = L As • (-VT) + ]TU B -F B 



(94) 



where the Onsager matrix (|93[) is positive definite, i.e. 



1 

T 



^ 



(95) 



-As i-AB 



The vectorial linear phenomenological relations ([Mf ac- 
count for the irreversibility due to the transport. 

Likewise, the pseudo- vectorial intrinsic rotation rate 
densities fl\ are related to the pseudo-vectorial torques 
T B through pseudo-vectorial linear phenomenological re- 
lations, i.e. 



/J L AB 



B J 



where the Onsager matrix has to satisfy, 



T { L ^4 



^ 



(96) 



(97) 



in order for the pseudo-vectorial quadratic form in the re- 
lation ([55)) to be positive definite. The pseudo-vectorial 
linear phenomenological relations (|96p are expressed ex- 
plicitly in terms of the intrinsic rotation rate densities fl\ 
and torques Ti as, 



'^a =J2{ L ab ■ (Pb x E) + L*g ■ (m B x B) 
= • (pb x E) + L™ M • (m B x B) 



where the Onsager matrix (|9T[) is positive definite, i.e. 



(98) 



1 
T 



I mm I mM 

AB AB 

I Mm I MM 

AB '-AB 



^ 



(99) 



The pseudo-vectorial linear phenomenological rela- 
tions (jMf account for the irreversibility due to the re- 
laxation. They also represent the classical counterpart of 
the spin-orbit coupling, since they couple the rotation rate 
density f2™ of the magnetisation with the rotation rate 
density of the matter J?™ . 

It is worth emphasising that although the scalar ([90)) . 
the vectorial (|M)) and the pseudo-vectorial (|98|) linear phe- 
nomenological relations are structurally independent, they 
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are coupled through the time evolution equations. The 
chemical affinities A a defined in relation (j74|) couple the 
scalar (|M|) and vectorial linear phenomenological re- 
lations. This coupling has interesting applications for spin- 
tronics, some of which were investigated in |13j . 

In the particular case of a continuous medium made 
of a single chemical substance A, the scalar linear rela- 
tions (j9"0)) vanish, the vectorial linear relations be- 
come, 



l ss -{-VT) + L sA -F A 
--l As -(-VT) + L AA -F A 



and the pscudo-vectorial linear relations (|98|) reduce to, 



njf = L^ 1 ■ ( PA x E) 4 
(pa x E) 



DM _ I Mn 
J ^A — L AA 



• (m A x B) 
h L^m . {mA x B ) 



4.2 Physical applications 

4.2.1 Lehmann effect 

In 1900, Lehmann [TJ] observed that droplets of liquid 
crystals, that have a chiral cholesteric structure, were set 
in rotation by a temperature gradient. The explanation 
given by Leslie |15j is now questioned |16j by recent ob- 
servations. Here, we show how our formalism accounts for 
the Lehmann effect. 

It is worth mentioning that the irreversible thermody- 
namics of a continuum of ncmatic liquid crystals has been 
examined by Miiller [5] and Muschik [?,?], even though 
these authors did not explicitly deduce the Lehmann ef- 
fect from their formalism. The "director" vector field they 
introduced corresponds to the preferred axis n of the con- 
tinuum of liquid crystals. In contrast to them, we do not 
enlarge the spatial manifold to account for n. Instead, we 
require the net electric polarisation P, which is a state 
field, to be collinear to n. 

We consider a uniform continuum made of identical liq- 
uid crystals that are oriented in the same direction. The 
continuum of liquid crystals is gyrotropic since the chiral- 
ity of the liquid crystals defines a preferred axis of unit vec- 
tor n. Each liquid crystal is helicoidal and made of iden- 
tical elements A that are dielectrics, i.e. q A = 0, and have 
electric dipoles p A orthogonal to the helix axis [17]. Thus, 
in the absence of an external interaction, the liquid crys- 
tals have no net electric polarisation. On a macroscopic 
scale, the continuum is homogeneous, i.e. V fi A = 0. We 
assume that the viscosity can be neglected, i.e. Vv^ = 
and V <jJ A = 0. The liquid crystals are trapped in an ex- 
perimental set-up such that they have no translational 
motion, i.e. j A = 0. Moreover, a temperature gradient 
V T is applied along a direction that is different from the 
preferred axis of the crystals. 

In such a case, the linear phenomenological rela- 
tion (|100p reduces to, 



where the temperature gradient V T induces a electric 
field E and the electric dipoles p A rotate in an asymmet- 
ric manner in order to lower the Debye energy — p A • E. 
Thus, a net electric polarisation P is generated along the 
preferred axis n, i.e. 

P = n A pa = n A (pa ■ n) n , (103) 

and the linear relation (|102|) is recast as, 

PVE= (n A L^ A -L As ) ■ VT . (104) 

The spatial symmetry requires the electric polarisation 
force P V E and the thermal force — V T to be collinear, 
i.e. 

P V E = n A L AA L As V T , (105) 

(101) where l A \ = L~ A \ \ and l As = L As 1. Thus, the phe- 
nomenological relation (| 105[) can be recast as, 

PVE = Ar i/ ifc B Vr, (106) 

where A is a dimcnsionlcss parameter given by, 



(100) 



A 



^aa Las 



(107) 



In the absence of an applied magnetic field, i.e. B = 0, 
the electric field E is expressed in terms of the electric 
potential V as, 

E = -W. (108) 

The electric polarisation force density P V E satisfies the 
vectorial identity, 



P V E = (P • V) E + P x (V x E) = (P • V) E 



p ■ 



(P0E) + (-V-P)E 



(109) 



where we used Faraday's law (]108[) and the index p de- 
notes that there is a dot product between the covariant 
differential operator V and the electric polarisation P. 
The term Vp • (P E) in the identity (]109[) corresponds 
to a surface contribution after integration over the volume 
of the liquid crystal continuum [18j . Thus, it can be ne- 
glected in the bulk of the continuum where the electric 
polarisation force density P V E is expressed as, 



P VE = (- V -P)E 



(110) 



Moreover, using the definition of the bound electric charge 
density gp, i.e. 

g p = -V-P, (111) 

the electric polarisation force density (|110p is recast 
as [TBI, 

PVE = g P E, (112) 

which shows that it is the analog of the electric part of the 
Lorentz force density for bound charges. Substituting the 
expression (|112p in the linear relation (|106[) . the latter can 
be expressed as a Seebeck effect for bound electric charges 
in an electric insulator, i.e. 



Pa VE 



-aa 



VT 



(102) 



E = £ P VT 



(113) 
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where the Seebeck coefficient for bound electric charges, 

Xua ks 

£ P = . 

<7P 

Moreover, the liquid crystals are subjected to an external 
electric torque ([57)1 that sets them in rotational motion. 
Thus, the external torque density is given by, 

r cxt = riA T ™ = P x E . (114) 

Using the equation (|113l) for the Seebeck effect, the torque 
density r ext is recast as, 

r cxt = £p (P x VT) . (115) 

The rotation occurs in the plane spanned by the electric 
polarisation P and the temperature gradient V T around 
the local centre of mass of the liquid crystals. Taking into 
account the absence of a matter current, i.e. = 0, and 
substituting the expression (|115p for the torque density 
into Newton's seond law in intrinsic rotation (|17p , the lat- 
ter accounts for the Lehmann effect, i.e. 

u, = ^(PxVT), (116) 

where / is the intrinsic angular mass density of the liquid 
crystals. The intensity of this effect is expressed as 

9 = £pPVT sm 6) ; (H7) 

where 9 is the angle between the electric polarisation 
P and the temperature gradient VT, i.e. P x VT| = 
P V T sin 9. 



4.2.2 Electric Lehmann effect 

In 1974, de Gennes [19] noticed that liquid crystals were 
set in rotation by an electric field. This effect is called 
the electric Lehmann effect [50] since it corresponds to 
a Lehmann effect where the driving force is an electric 
force instead of a thermal force. Here, we show how our 
formalism accounts for the electric Lehmann effect. 

We consider the same continuum of liquid crystals as in 
Sec. 14.2.11 In contrast to the Lehmann effect, the electric 
Lehmann effect is an isothermal effect. The temperature 
gradient V T is replaced by an electric potential gradient 
V V that is applied along a direction that is different from 
the preferred axis of the liquid crystals. By analogy with 
the Lehmann effect, the electric field must be inhomoge- 
neous in order to generate an electric polarisation force 
density (fTT2"]l . 

In the presence of an applied electric potential gradi- 
ent V V, the electric dipoles rotate in an asymmetric 
manner in order to lower the Debye energy p^- V V, which 
yields a net electric polarisation P along the preferred axis 
n given in equation (|103|) . Moreover, the liquid crystals 
are subjected to an external electric torque density (|114[) , 
that sets them in rotational motion. In the absence of a 



magnetic induction field B, Faraday's law (|108p implies 
that the torque density (|114|) yields, 

r ext = -PxVU. (118) 

The rotation occurs in the plane spanned by the electric 
polarisation P and the electric potential gradient V V 
around the local centre of mass of the liquid crystals. 
Taking into account the absence of a matter current, i.e. 
= 0, and substituting the expression (|118p for the 
torque density into Newton's second law in intrinsic ro- 
tation (|T7|) . the latter accounts for the electric Lehmann 
effect, i.e. 

u = -j(PxVV) . (119) 
The intensity of this effect is expressed as 

PW 

9 = — sin (2 , (120) 

where 9 is the angle between the electric polarisation P 
and the electric potential gradient V V, i.e. |P x V V \ = 
P\7V sin9. 

It is useful to point out that Quincke [21] observed in 
1896 that small dielectric spheres in suspension in a liquid 
placed inside a parallel-plate capacitor are set in rotational 
motion by charging the capacitor. The intrinsic rotational 
dynamics of the Quincke effect is described by the same 
time evolution equations (|119p and (|120p as the electric 
Lehmann effect. 



4.2.3 Relaxation of electric dipoles 

The first model of the relaxation of electric dipoles pa 
is due to Debye [22]. Here, we show how our formalism 
accounts for this relaxation. 

We consider a homogeneous sample made of a single 
chemical substance A, consisting of electric dipoles pa in 
the absence of a magnetic induction field, i.e. B = 0. In 
the frame of sample, the chemical substance A has no 
translational motion, i.e. ]a = 0. 

The absence of a magnetic induction field, i.e. B = 0, 
implies that the linear phenomenological relation (jlOip 
reduces to, 

f2% = LTX ■ (PA x E) . (121) 

The process is irreversible and dissipativc, which means 
that the Debye energy — p^ • E has to diminish. This 
implies that the intrinsic rotation of the local element, 
to which the electric dipolc p^ is attached, occurs in 
the plane spanned by p^ and E. Thus, the rotation 
rate density 12™ is collinear to the torque pa X E and 
L™™ = LTa 1 - Substituting the relation (fHTj) into the 
time evolution equation Q33p for p^ yields an equation ac- 
counting for the Debye relaxation of the electric dipoles 
Pa in the presence of an electric field E, i.e. 

p A = -a A PAX (Pa x E) , (122) 

where 

a A = n-/ L™ , 
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is a phcnomcnological friction coefficient. 

The electric field E is an effective field that is defined 
with respect to the local infinitesimal system. It consists 
of two contributions, i.e. 



E = E ( 



E 1 



(123) 



The first contribution E oxt is due to an external field ap- 
plied on the whole system. The second contribution E lnt 
is due to the dipolar interaction with the infinitesimal sys- 
tems that are in the neighbourhood of the local system. 
Note that this contribution is internal to the whole system, 
but external to the local infinitesimal system. In the neigh- 
bourhood of the local system, the Debye energy density 
due to the dipolar interaction between the local systems is 
proportional to the magnitude of the spatial variation of 
the electric polarisation. Thus, the Debye energy density 
can be expressed as, 



P E 



P E c 



D P VP 



(124) 



where D is a phcnomcnological parameter. Note that the 
spatial variation is expressed in terms of a Laplacian V 2 
and not a gradient V, since otherwise the contributions 
of local neighbouring systems located on opposite sides of 
the local system would cancel out. The expressions (|123D 
and (|124|) imply that the electric field E mt is given by, 



E 



hit 



D V 2 P 



For a single homogeneous sample made of a single chemi- 
cal substance A, consisting of electric dipoles pa, the elec- 
tric polarisation P = uaPa- Thus, the time evolution 
equation (|122p becomes, 

p A = -a A PA x (p A xE oxt ) (126) 
- D a A Pa x (p A x V 2 (n A Pa) ) , 

where the terms on the RHS describe the relaxation due 
to an external electric field E oxt and to the interaction 
with the electric dipoles in the neighbourhood of the local 
system. 



4.2.4 Relaxation of magnetic dipoles and magnetisation 
waves 

The first model of relaxation of magnetic dipoles in 
the presence of a magnetic induction field B is due to 
Landau and Lifshitz [22- Later, Gilbert J2U developed 
a Lagrangian formulation of this damping. Here, we show 
how our formalism accounts for this relaxation and for the 
magnetisation waves, which are the classical counterpart 
of the "spin waves" [25] . 

We consider a homogeneous sample made of a single 
chemical substance A, consisting of magnetic dipoles 
in the absence of an electric field, i.e. E = 0. In the frame 
of sample, the chemical substance A has no translational 
motion, i.e. ]a = 0. 



The absence of an electric field, i.e. E = 0, implies that 
the linear phenomenological relation (jlOip reduces to, 



f2 



M 



MM 

-AA 



(m A x B) 



(127) 



The process is irreversible and dissipativc, which means 
that the Larmor energy — • B has to diminish. This 
implies that the intrinsic rotation occurs in the plane 
spanned by m A and B. Thus, the rotation rate density 
ilj/ 1 is collincar to the torque x B and = 
Laa t. Substituting the relation (fT2T|) into the time evo- 
lution equation (|40[) for m A yields the Landau-Lifshitz 
equation |24| accounting for the relaxation of the mag- 
netic dipoles m A in the presence of a magnetic induction 
field B, i.e. 



ih A = 1a mA x B - [3 A m A x (m A x B) 



(128) 



where 



-1 t MM 

l A ^AA 



is a phenomenological friction coefficient. Note that the 
relaxation of the magnetic dipoles is due to a mix of bod- 
ily rotation [26j composed with Neel-type magnetic relax- 
ation [27] , 

The magnetic induction field B is an effective field that 
is defined with respect to the local infinitesimal system. It 
consists of two contributions, i.e. 



B = B' 



B 1 



(129) 



(125) The first contribution B oxt is due to an external field ap- 



plied on the whole system. The second contribution B 1 is 
due to the ferromagnetic interaction with the infinitesimal 
systems that are in the neighbourhood of the local system. 
In the neighbourhood of the local system, the Larmor en- 
ergy density due to the ferromagnetic interaction between 
the local systems is proportional to the magnitude of the 
spatial variation of the magnetisation [28j . Thus, the Lar- 
mor energy density can be expressed as, 



M B = M B AM V M 



(130) 



where A is a phenomenological parameter called the stiff- 
ness constant [35] . The expressions (|129p and (|130[) imply 
that the magnetic induction field B mt is given by, 

AV 2 M . 



B 



int 



(131) 



For a single homogeneous sample made of a single chem- 
ical substance A, consisting of magnetic dipoles m A , the 
magnetisation M = n A va. A . Thus, the Landau-Lifshitz 
equation (|128[) becomes, 

rh A = 1A m A x B ext - (3 A m A x (m A x B cxt ) 

+ A 1A m A x V 2 (n A m A ) (132) 

- APavha x (m A x V 2 (n A m A )^ ■ 

The first and second terms on the RHS of the Landau- 
Lifschitz equation (|132p describe respectively the preces- 
sion and relaxation due to an external magnetic induc- 
tion field B cxt . The third term describes magnetisation 
waves [55] and the fourth term describes the relaxation of 
magnetisation waves. 
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4.2.5 Thermally driven magnetisation current 

The relaxation of magnetic dipoles in a ferroelectric metal, 
established in Sec. 14.2.41 rests on the assumption that 
the chemical units carrying the magnetic dipoles are at 
rest. Here, we consider a ferromagnetic conductor made 
of a fixed lattice and conduction electrons. The magnetic 
dipoles are carried by the core electrons of the lattice and 
the conduction electrons (e.g. 'd-f and 's-p' electrons re- 
spectively). On the scale of interest, the core electrons 
of the lattice are described as a continuum of chemical 
type A at rest and the conduction electrons as a fluid 
of chemical type B. The core electrons carry a magnetic 
dipole xua and the conduction electrons a magnetic dipole 
m B . A quasi- uniform and constant magnetic induction 
field B and temperature gradient V T arc applied on the 
ferromagnetic conductor. No electric field is applied, i.e. 
E = VV = 0. 

Since the core electrons are at rest on the scale of 
interest, the current density of the core electron contin- 
uum vanishes, i.e. j A = 0. In our formalism, we neglect 
the "chemical reactions" a between the core and conduc- 
tion electrons, i.e. u) a = 0. Thus, the time evolution equa- 
tion (|4U|) for the magnetic dipoles of the core electron 
continuum A reduces to, 



n A rh A = j A n A (m A x B) - m A x fl™ 
+ Jab n A (m A x n B m B ) . 



(133) 



The absence of electric field, i.e. E = 0, implies that the 
linear phenomenological relation ([55)1 for the core electron 
continuum reduces to, 



= • (m A x B) + • (m B x B) . (134) 

The effect on the rotation of the core electrons of the mag- 
netic torque niflxBof the conduction electrons is negligi- 
ble compared to the effect of the magnetic torque x B 
of the core electrons, i.e. L A b -C \-aa, which yields the 
phenomenological equation (|127p . The magnetic induction 
field B consists of two contributions as established in equa- 
tion (|129[) where the magnetisation M = n A m A + n B m B , 
i.e. 



B = B CIt +iV 2 (n A m A + n B m B ) 



(135) 



Following the same lines of thought as in Sec. 14.2.41 the lo- 
cal time evolution equation (| 1 33f> for the magnetic dipoles 
of the core electron continuum is found to be, 

rh A = 7a (m A x B oxt ) - f3 A m A x (m A x B oxt ) 
+ Aj A m A x V 2 (n A m A + n B m B ) 

- A p A m A x (mA x V 2 (n A m A + n B m B ) 
+ j AB (m A x n B m B ) ■ (136) 



The last term on the RHS of the equation (| 136[) is ex- 
pected to lead to the magnetisation transfer torque, which 
is the classical counterpart of the "spin transfer torque" , 
generated by the magnetic dipoles m B of the conduction 



electron fluid on the magnetic dipoles of the core elec- 
tron continuum [30]. Note that the time evolution of the 
magnetic dipoles m A of the core electron continuum has 
no explicit dependence on the temperature gradient be- 
cause these electrons do not undergo transport. However, 
as shown below, the dynamics of the magnetic dipoles m B 
depends on the temperature gradient. 

Since we neglect the "chemical reactions" a between 
the core and conduction electrons, i.e. u) a = 0, the time 
evolution equation (|40p for the magnetic dipoles mg of 
the core electron continuum B is given by, 

n B m B = 7s n B (m B x B) - m B x - (j B • V) m B 
+ 1b A n B (m B x n A m A ) . (137) 

Following the same procedure as for the core electrons, 
we recast the local time evolution equation (|133p for the 
magnetic dipoles of the conduction electron fluid as, 

m B = 7s (m B x B cxt ) - f3 B m B x (m B x B cxt ) 

+ A-f B m B xV 2 (n A m A + n B m B ) (138) 

- A/3 B m B x (m B x V 2 {n A m A + n B m B ) 

- (js ' V) m B + j BA (m B x n A m A ) . 

In order to find an explicit expression for the third term on 
the RHS of the time evolution equation (| 138(1 , we describe 
the transport of the conduction electron fluid in the pres- 
ence of a uniform magnetic induction field B and a tem- 
perature gradient V T. In the absence of an electric field, 
i.e. VF = 0, the linear phenomenological relation (jlOOj) 
reduces to, 

Jb = - — <r B -£ B -VT, (139) 
Qb 

where the isothermal electric conductivity tensor <r B and 
the thcrmopowcr tensor e B are defined respectively as [5] , 



0" b = <Z B L BB , 
1 

IB 



£b = — \- BB ■ L Bs 



(140) 



and satisfy the following symmetries [3~T] . 

- B (B)= CT |(-B) , 
e B (B) = 4(-B) , 

where the exponent T stands for transpose. Thus, these 
phenomenological tensors can be split into symmetric and 
antisymmetric parts according to, 



a B ■ a = <7 B || a+ <r B ± (B x a) , 



e B • a = e B ii a + e B ± B x a 



(141) 



where the magnetic induction field B = B 6 and a is an 
arbitrary vector. By symmetry the collinear factors <7 B || 
and e B || are independent of the intensity B of the mag- 
netic induction field. The Hall and Nernst effects [9ll3Tj 
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require the orthogonal factors <t B j_ and e B j_ respectively 
to be inversely proportional to B, i.e. cr B ±_ oc B' 1 and 
e B ± oc B -1 . Using the splittings (|14ip . the linear phe- 
nomcnological relation (|139[) becomes, 



Jb 



a B\\ £ B\\ 

Qb 



vr 



a B\\ £ B± + &B± £b\\ 
IB 



B x VT 



(142) 

where the first term on the RHS describes the Soret ef- 
fect [5] for the component of the magnetic induction field 
B that is collinear to the temperature gradient V T and 
does not affect the transport, and the second term de- 
scribes the Ettingshausen effect [5] per unit of electric 
charge qs for the components of the magnetic induction 
field B that are orthogonal to the temperature gradient 
VT and affect the transport. 

Since the core electrons A do not participate to the 
transport, the magnetisation current density tensor jM is 
entirely due to the transport of the conduction electrons 
B. Using the relations ((37)) and (|142[) . jM is found to be, 



JM = m B V T 

qB 

vb\\ £b±. + &BX £b\\ 



Qb 



m B (B x VT) 



(143) 



Substituting the linear phenomenological equa- 
tion (|142[) into the local time evolution equation (|138[) for 
the magnetic dipoles ms of the conduction electron fluid 
B, the latter becomes, 

m B = 7s (m B x B oxt ) - (3 B m B x (m B x B cxt ) 

+ A^f B m B xV 2 (n A m A + n B m B ) (144) 
- A/3b m B x (m B x V 2 (n A m A + n B m B ) 



®B II & B II 

+ 7 B a (m B x n A m A ) H - V T ■ V m B 

qB n B 

+ *B 1 e B± + a B± e Bl /^^.^^ 
q B n B \ ) 

where the first and second terms on the RHS account re- 
spectively for the precession and the relaxation of the mag- 
netic dipoles m B . The third and fourth terms account re- 
spectively for the magnetisation waves and the relaxation 
of the magnetisation waves. The fifth term accounts ac- 
counts for the interaction between the magnetic dipoles 
and m B of the core and conduction electrons respec- 
tively. The sixth and seventh terms account for the trans- 
port of m B and describe the magnetisation accumulation, 
which is the classical counterpart of the "spin accumula- 
tion" , generated by the temperature gradient V T in the 
presence of a magnetic induction field B. The sixth term 
describes the magnetisation accumulation collinear to the 
temperature gradient and the seventh term describes the 
magnetisation accumulation orthogonal to the tempera- 
ture gradient and the magnetic induction field. 

The times scales associated to the precession and the 
relaxation of the conduction electrons are much smaller 



than the time scales associated to the magnetisation ac- 
cumulation and the magnetisation transfer torque. Thus, 
on the latter time scales, the first, second, third and fourth 
terms on the RHS of the time evolution equation (|144p can 
be neglected, i.e. 

m B = j BA (m B x n A m A ) + B|1 £g|1 V T • V m B 



£B±_ + &B1. SB\\ 

qB n B 



qBn B 
B x V T ) • V m B 



(145) 



4.2.6 Electrically driven magnetisation current 

The time evolution of the magnetic dipoles of a ferro- 
magnetic conductor, made of a fixed lattice and conduc- 
tion electrons, in the presence of a quasi-uniform magnetic 
field B and a temperature gradient V T was established 
in Sec. 14.2.51 Here we consider an isothermal conductor, 
i.e. VT = 0, in the presence of a quasi-uniform electric 
field, i.e. E = — V V. On the scale of interest, the core 
electrons of the lattice are described as a continuum of 
chemical type A at rest and the conduction electrons as a 
fluid of chemical type B. The core electrons carry a mag- 
netic dipole and the conduction electrons a magnetic 
dipole m B . 

The core electron continuum is at rest and the local 
time evolution equation (| 136[) for the magnetic dipoles 
m A is independent of the transport. The time evolution 
for the magnetic dipoles m B of the conduction electron 
fluid is given by the equation (|138|) . 

In order to find an explicit expression for the third 



term on the RHS of the time evolution equation (|138[) , we 
describe the transport of the conduction electron fluid in 
the presence of a magnetic induction field B and a electric 
potential gradient V V. In the absence of a temperature 
gradient, i.e. VT = 0, the linear phenomenological rela- 
tion (|100p reduces to, 



j B = - — er B • V V 

qB 



(146) 



Using the first splitting (|14ip , the linear phenomenological 
relation (|146p becomes, 



P£|| 

qs 



vu- 



qB 



B x V V 



(147) 



where the first term on the RHS is Ohm's law [5] per unit 
of electric charge q B for the component of the magnetic 
induction field B that is collinear to the electric potential 
gradient V V and does not affect the transport, and the 
second term describes the Hall effect [S] per unit of electric 
charge q B for the components of the magnetic induction 
field B that are orthogonal to the electric potential gradi- 
ent V V and affect the transport. 

Since the core electrons A do not participate to the 
transport, the magnetisation current density tensor jM is 
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entirely due to the transport of the conduction electrons 
B. Using the relations ((37|) and (|147|) . jM is found to be, 

Jm = - — m B &VV- — m B 0fB x V v) . (148) 

Substituting the linear phenomcnological equa- 
tion (|147[) into the local time evolution equation (|138[) for 
the magnetic dipoles ms of the conduction electron fluid 
B, the latter becomes, 

m B = 7s (m B x B cxt ) - /3 B m B x (m B x B cxt ) 

+ A"f B m B x V 2 {n A m A + n B m B ) (149) 
- A fl B m B x (m B x V 2 (n A m A + n B m B ) 

°BII 



7 B A (m B x n A m A 

<JB 



q_B n B 



qB n B 
BxW) • Vm B 



■ VV ■ Vm B 



where the sixth and seventh terms account for the trans- 
port of m B and describe the magnetisation accumulation 
generated by the electric potential gradient V V in the 
presence of a magnetic induction field B. The seventh term 
describes the magnetisation accumulation collinear to the 
electric potential gradient and the fifth term describes the 
magnetisation accumulation orthogonal to the electric po- 
tential gradient and the magnetic induction field. 

On the time scales associated to the magnetisation 
accumulation and the magnetisation transfer torque, the 
time evolution equation (|149[) reduces to, i.e. 



me = 7ba (m B x n A m A ) 

&B. 



a B\\ 



qB n B 



qB n B 

B x VV) ■ Vm B 



VV ■ Vm B 



(150) 



4.2.7 Thermodynamics of magnetic vortices 

The magnetic counterpart of electrically polarised liq- 
uid crystals are magnetic vortices or skyrmions [52"] . 
Skyrmions were observed recently |33| in the insulating 
fcrromagnet Cu20Se03. Our formalism predicts the pre- 
cession and relaxation of magnetic vortices in the presence 
of a temperature gradient. 

We consider magnetic vortices that can be treated as 
a continuum. They form a gyromagnctic medium with a 
vorticity axis of unit vector n. We assume that n is uni- 
form. The magnetic vortices are made of ions and core 
electrons and are considered as an electrically neutral sub- 
stance of type A, i.e. q A = 0. In the absence of an exter- 
nal interaction, the vortices have no net magnetisation. 
On a macroscopic scale, the continuum is homogeneous, 
i.e. V /j, a = 0. We assume that the viscosity can be ne- 
glected, i.e. i.e. Vva = and V u> A = 0. Moreover, since 
the magnetic dipoles are carried by the electrons, that 
have no intrinsic angular mass, we do not consider the 
rotation of the matter but we take into account only the 
rotation of the magnetisation. In an electrical insulator, 



the magnetic vortices have no translational motion, i.e. 
j A = 0. Furthermore, a temperature gradient V T is ap- 
plied along a direction that is different from the vorticity 
axis of the magnetic vortices. 

In this case, the linear phenomenological relation (|100[) 
reduces to, 



m A VB= (L- A \-L As ) -VT , 



(151) 



where the temperature gradient V T induces a magnetic 
induction field B and the magnetic dipoles rotate in 
an asymmetric manner in order to lower the Larmor en- 
ergy — m A ■ B. Thus, a net magnetisation M is generated 
along the vorticity axis n, i.e. 

M = n A m A = n A (m A • n) n , (152) 

and the linear relation (|151[) is recast as, 

MVB=(n A L A 1 A -L As ) -VT . (153) 

The spatial symmetry requires the magnetisation force 
M V B and the thermal force — V T to be collinear, i.e. 



M V B = n A L- A \ L As V T , 



(154) 



where \- AA = L A l A 1 and L^ s = L As 1. Thus, the phe- 
nomenological relation f|l 54[) can be recast as, 



MVB = Xn A k B VT , 



(155) 



where A is a dimcnsionlcss parameter by the relation (|107[) . 
The magnetisation force density MVB satisfies the vec- 
torial identity, 

MVB = V(M-B)-BVM-Bx(VxM) (156) 
= — Vm • (M Q B — (M • B) 1) + ( V x M) x B 



where we used Thomson's law, i.e. 



V • B = , 



(157) 



and the index m denotes that there is a dot product be- 
tween the covariant differential operator V and the mag- 
netisation M. The term V M • (M B - (M • B) 1) in the 
identity (|156|) corresponds to a surface contribution after 
integration over the volume of the continuum of magnetic 
vortices [TB]. Thus, it can be neglected in the bulk of the 
continuum where the magnetisation force density MVB 
is expressed as, 



MVB=(VxM)xB 



(158) 



Moreover, using the definition of the bound electric cur- 
rent density jm, i-c 

jm = V x M , (159) 
the magnetisation force density (|158[) is recast as [18], 

M V B = jm x B , (160) 
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which shows that it is the analog of the magnetic part of 
the Lorentz force density for bound currents. Substituting 
the expression (|160[) in the linear relation (|1 54[) , the latter 
becomes, 

j M x B = \n A k B VT , (161) 
which implies that 

Jm x (j M x B) = A n A k B (jm x V T) , (162) 

and that the thermally induced magnetic induction field 
B is orthogonal to the bound current jmj i-e. 

Jm • B = . 

Using the vectorial identity 

Jm x (j M xB) = (j M • B) j M - Jm B — _ Jm B i 

the linear relation (|162p describing the magnetic field B 
induced by a temperature gradient V T on magnetic vor- 
tices jM can be recast as a magnetic Seebeck effect for 
bound currents in an electric insulator, i.e. 

B = £ M xVT, (163) 

where the magnetic Seebeck vector for bound electric cur- 
rents, 

£m = - Att-a k B j M ■ 

Note that the magnetic Seebeck effect (|163[) is the mag- 
netic analog of the Seebeck effect (|113p in an electric in- 
sulator. 

Substituting the expression (|163|) for the induced mag- 
netic field B into the local time evolution equation of the 
magnetic dipoles (|128|) . we obtain an expression for the 
thermally induced dynamics of the magnetic vortices, i.e. 

= 7^ m A x (e M x V T) (164) 

- /3 A m A x (m A x (e M xVT)j , 

where the first term and second terms on the RHS de- 
scribe respectively the precession and the relaxation of 
the magnetic dipoles of the core electrons of the magnetic 
vortices due to magnetic field B induced by a temperature 
gradient V T. 

4.2.8 Thermally driven magnetisation waves 

We consider a ferromagnetic insulator made of a fixed lat- 
tice with core electrons that carry a magnetic dipole m A . 
A constant external magnetic induction field B cxt and a 
temperature gradient V T are applied. The time evolution 
of the magnetic dipoles is given by the Landau-Lifschitz 
equation (|128[) . 

The magnetic induction field B appearing in the time 
evolution equation (|128[) is an effective field that is defined 
with respect to the local infinitesimal system. It consists 
of three contributions, i.e. 

g = gcxt +B int +B ind _ ( 165 ) 



The first contribution B oxt is due to an external field ap- 
plied on the whole system. The second contribution B mt 
is due to the ferromagnetic interaction with the infinites- 
imal systems that are in the neighbourhood of the local 
system. For a core electron continuum, the magnetisation 
M = n A m A and the expression for B lnt follows from 
equation (|131[> . i.e. 

B int = AV 2 {n A m A ) . (166) 

The third contribution B lnd is induced by the temper- 
ature gradient V T. For a core electron continuum, the 
expression for B lnd follows from equation (|163p . i.e. 

B ind = e M xVT, (167) 

where the bound current density yields, 

Jm = V x (^m^) . (168) 

Note that the magnetic induction field B lnd is induced 
only in the presence of a bound current density jM- This 
current density arises in the presence of magnetisation 
waves generated by the applied external magnetic field 

B oxt 

Substituting the expressions (|166p and (|167p for the 
different contributions to the effective magnetic field B 
into the Landau-Lifschitz equation (|128p . the latter be- 
comes, 

rh A = j A (m A x B oxt ) - f3 A m A x (m A x B cxt ) 

+ j A m A x (AV 2 (n A m A ) + e M x VT) (169) 

- P A m A x ^m A x (A V 2 (n A m A ) +£ M xVTj| . 

where the third and fourth terms describe respectively 
how magnetisation waves and their relaxation are driven 
by a temperature gradient V T. The effect is maximal 
when the temperature gradient is orthogonal to the bound 
current jM generated by the applied external magnetic 
field B oxt . 



5 Conclusion 

The thermodynamics of irreversible processes is con- 
sidered for an electrically charged continuous medium 
containing spontaneous electric and magnetic dipoles in 
the presence of electromagnetic fields. Expressing the 
extensive matter state fields in terms of their chemi- 
cal constituents yields explicit expressions for the cur- 
rent densities. Three types of irreversible terms are de- 
rived from these expressions. These are scalar, vectorial 
and pseudo-vectorial terms that describe respectively ir- 
reversible chemical processes, irreversible transport pro- 
cesses and irreversible relaxation processes. These pro- 
cesses are coupled through the time evolution equations. 
Note that with such an approach, the mathematical struc- 
ture of the irreversible thermodynamics is uncovered phys- 
ically without imposing it a priori using irreducible repre- 
sentations of the Euclidean group. 
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As an illustration of our formalism, wc describe no- 
tably the Lehmann and electric Lehmann effects, the re- 
laxation of the electric and magnetic dipoles. In particular, 
we are able to predict the effect of a temperature gradient 
on the time evolution of the magnetic dipoles of conduc- 
tion electrons interacting with core electrons, which leads 
to precise expressions for the thermal and electronic mag- 
netisation accumulations. We also predict the precession 
and relaxation of magnetic vortices induced by a tempera- 
ture gradient in the absence of an applied magnetic induc- 
tion field, which is very innovative. Finally, in the presence 
of an applied magnetic induction field, we show explicitly 
how a temperature gradient drives magnetisation waves. 

The authors would like to thank Fran§ois A. Reuse for theo- 
retical guidance. 
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